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On a class of semilinear fractional elliptic equations 
involving outside Dirac data 


Huyuan ChenQ Hichem Hajaiejl Ying \Van 'll 


Abstract 

The purpose of this article is to give a complete study of the weak solutions of the fractional 
elliptic equation 

(— A) a u + u p = 0 in £?i(ejv), 

u = Sq in R N \ -Bi(eAr), 

where p > 0, (—A)“ with a £ (0,1) denotes the fractional Laplacian operator in the principle 
value sense, B\(e n) is the unit ball centered at ejss = (0, • • • ,0,1) in R w with N > 2 and do is 
the Dirac mass concentrated at the origin. We prove that problem (ED admits a unique weak 
solution when p > 1 + Moreover, if in addition p > AA|, the weak solution vanishes as 
a —> 1 “. We also show that problem m doesn’t have any weak solution when p £ [0,1 + ^]. 
These results are very surprising since there are in total contradiction with the classical setting, 
i.e. 

—Au + u p = 0 in i?i(ejv), 

u = Sq in \ -Bi(ejv), 

for which it has been proved that there are no solutions for p > . 
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1 Introduction 

Fractional PDEs have gained tremendous interest, not only from mathematicians but also from 
physicists and engineering, during the last years. This is essentially due to their widespread domains 
of applications. In fact the fractional Laplacian arises in many areas including medicine m, bio¬ 
engineering MM, 2] [22], relativistic physics [II 1171 ITS]. Modeling populations [29], flood flow, 
material viscoelastic theory, biology and earthquakes. It is also particularly relevant to study some 
situations, in which the fractional Laplacian is involved in PDEs, featuring irregular data such that 
those phenomena describing source terms which are concentrated at points. In our context, the 
source is placed outside the unit ball B i(eAr). This generates long-term interactions and short¬ 
term interactions, described by the nonlocal operator (—A)" and the nonlinear absorption u p 
respectively. (—A) Q has also a probabilistic interpretation, related to the above one. It is the 
a—stable subordinated infinitesimal killed Brownian motion. 

Let -E>i(e,y) be the unit ball in W N (N > 2) with center ejv = (0, • • • , 0,1) and So be the Dirac 
mass concentrated at the origin. Our main objective in this article is to investigate the existence, 
nonexistence and uniqueness of positive weak solutions of the semilinear fractional equation 

(— A) a u + u p = 0 in Bi(e]y), 

u = do in W N \ L?i(ejv), 

where p > 0 and the fractional Laplacian (—A)“ with a £ (0,1) is defined by 

(—A ) a u(x) = CN, a lim (—A )^u(x), 

e—>0+ 
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where 


-l 


CN,a 


with z = (z i,---, zjv) £ R N and 


1 — cos(zi) 
\z\ N+2a 



(-A)>(x) = - [ 

JR N \B e (x) 


u(z ) — tt(x) ; 
|._ x |7V+2a^ 


( 1 . 2 ) 


In 1991, a fundamental contribution to semilinear elliptic equations involving measures as 
boundary data is due to Gmira and Veron (15], where they studied the existence and uniqueness 
of weak solutions for 


—A u + h(u ) = 0 in SI, 
u = p on <912, 


(1.3) 


where 12 is a bounded C 2 domain and p is a bounded Radon measure defined in <912. A function u 
is said to be a weak solution of (11.31) if u £ L 1 (12) J h(u) £ L 1 (12, / 9<2x) and 


[ [«(-A)£ + h{u)£\dx = f ^^-dp(x), € C 0 L1 (12), (1.4) 

Jn Jon an x 

where p(x) = dist(x, <912) and n x denotes the unit inward normal vector at a point x. Gmira and 
Veron proved that the problem (11.31) admits a unique weak solution when h is a continuous and 
nondecreasing function satisfying 



[h(s) — h(—s)]s 1 


JV+l 

N - 1 ds < +oo. 


(1.5) 


The weak solution of (11.31) is approached by the classical solutions of (11.31) when p is replaced by 
a sequence of regular functions {p n }, which converge to p in the distribution sense. Furthermore, 
they showed that there is no weak solution of m when p = S xo with xq £ <912 and h(s) = |s| p 
with p > Later on, this subject has been vastly expanded in recent works, see the papers of 

Marcus and Veron [ 231122123122 ], Bidaut-Veron and Vivier [3] and references therein. 

In the fractional setting, the equivalent of (11.31) when p = 5 Xo has been considered in [?], where 
the authors proved that the weak solution of 

(-A ) a u + g(u) = k^P- 
u = 0 


in 12, 

in 12 c 


( 1 . 6 ) 


is approximated by the weak solutions, as t —» 0 + , of 


(-A) a u +g(u) = kt a 5 Xo+tfixo in 
u = 0 in 


12 , 

12 c . 


Qot 

More precisely, in the fractional setting, gJf 1 plays the same role of u = S xo on <912 in m- Our 
purpose in this article is to study the solution of 

(— A) a u + u p = 0 in ^i(ejv) 

when the exact Dirac mass concentrated at the origin is considered. Our main idea is to make use 
of nonlocal properties of the fractional Laplacian to move the Dirac mass at —te n when t —> 0 + 
and we then proceed by approximation techniques. Before giving our main results, we must first 
give an appropriate definition of weak solution of (11.11) . It is then worth to mention two important 
results. The equation 


(-A ) a u + uP = 0 
u = f 


2 


in B i(ejv), 
in R N \B 1 (e N ), 


(1.7) 








where / £ CoQKf^ \ B i(ctv)), admits a unique classical solution Uf, see [6, Theorem 2.5]. Further¬ 
more, let Uf = Uf in B i(ejv) and Uf = 0 in W N \ Bi(ei\r), then Uf is the unique classical solution 
of 

(~A) a u + u p = c N , a f R N\ Bl{eN) \ x Jy\N+-i a dy in B^{e N ), (1 g) 

u = 0 in R JV \Bi(e JV ) 


and satisfies the identity: 

/ [u(x)(-A) a £(x) + u p (x)C(x)] dx = c Nta [ [ 

J(eiv) J B-\ (eiv) 


Bi(ejv) JR^B^eiy) \ x ~ V\ 


Z{x)f{y) 

_ ni\N-\-2cX 


dydx , 


for any £ £ Cfi°(Bi(e]\r)). Let us mention that C^°(Bi(ei\r)) is the space of test functions £ £ 
C°°(R n ) with support in Bi(cn)- 

Inspired by above identity, we give the definition of weak solution to (11.11) as follows. 

Definition 1.1 ITe say that u is a weak solution of 11.11 ) if u £ L l (B\(eN)), u p £ L\ oc (B\(e n)) 
and 


L 


Bi(ejv) 


[u(x)(— A)“£(x) + u p (x)f(x)] dx = 


' Bi (ejv) 


£(x)r 0 (x)dx, V££C' 0 oo (S 1 ( ejv )), (1.9) 


where 

r oW = ^L Vx£M Af \{0}. (1.10) 

It is well known that the definition of the weak solution heavily depends on the test functions 
space and the best function space is the one which enables us to get the ” strongest” weak solution, 
la fflj [lQ]j semilinear fractional equations with measures has been studied via the test functions 
space C C(M. N ) for a C 2 bounded open domain Q, where X ai n is the space of functions £ 

satisfying: 

(1) supp(£) C Q; 

(2) (—A) Q £(x) exists for all x e Q and |(—A)"£(x)| < C for some C > 0; 

(3) there exist p £ L 1 {£l,p a dx) and eo > 0 such that |(—A)“£| <p a.e. in D, for all e £ (0, eo], 

where p(x) = dist(x,dil). 

The test functions space Cfi°(B i(ejv)) has stronger topology than X QjBl ( ejv ) does, the weak solution 
in Definition 1.1 with test functions space 'A 0 .Bi(e N ) would be stronger than the one with test 
functions space Co°(Bi(ei\r)). It is then worth to mention that the test functions space C , Q°(l?i(eAr)) 
could not be replaced to the test functions space X aBl ( ejv ) in our setting. For example, £o := 
G q [1] £ X Q;Sl(ejv) , but ()1.9D does not hold for £o, where G a denotes the Green kernel of (—A) Q in 
-Bi(ejv) x B\ (e^v) and G a is the Green operator defined as 


G a [f](x) = [ G a (x,y)f{y)dy , / £ L 1 (Bi(e N ) , p a dx) . (1.11) 

J Bi(ejv) 


Let us state our existence result. 


Theorem 1.1 Assume that a £ (0,1) and p > 1 + jf. Then there exists a unique nonnegative 
weak solution u aiP of 11.11 ) such that for some c\ > 1, we have 

N+2a 

0 < u a , P (x) < ci\x\ p , Vx £ -Bi(ejv) (1.12) 

and 

| jV+2q N+2a 

— t p < Ua^tejy) < C\t P , Vt £ (0,1). (1.13) 
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Remark 1.1 (i) The existence results are very surprising as they are in total different from, the 
Laplacian case, where hi. 3 1) with p = 5 q has a weak solution only when p < • 

(ii) From \1.13\) . the singularity is only near the origin. We also notice that 

u ap(te at ) > —t~^ N+2a \ \/t £ ( 0 , 1 ), 

Cl 

which implies that the absorption nonlinearity u p plays a primary role in the equation cop. While 
the absorption nonlinearity always plays a second role in a measure framework. 

(Hi) The uniqueness cannot directly follow Kato’s inequality m Proposition 2.f] since it has 
been built in the framework of the test functions space X QiSl ( ejv ). In this paper, as mentioned 
above, Cff(B\(eN )) is the appropriate test functions space. This will give birth to a lot of technical 
difficulties to prove the existence, nonexistence and uniqueness of weak solutions of mn- 

ii p > 1 + ^, the weak solution m Q iP of m is approximated by the unique solution u s 
(s £ (0,1)) of 

(— A) a u + u p = 0 in B\(ei\f), 

AT (1.14) 

u = 5- seN in R N \B 1 (e N ). 

When p £ [0,1 + ^], we will prove that {u s } blows up everywhere in B\(cn) as s —> 0 + , 
therefore, we can deduce the nonexistence of weak solutions of CEO when p < 1 + ^. More 
precisely, we have the following results. 

Theorem 1.2 Assume that a £ (0,1) and 0 < p < 1 + ^. Then problem 11.11) does not have any 
weak solution. 


In the proof of Theorem 11.21 we will first need to prove the crucial estimate 

<p(*) >c 2 |x|-( jv+2q ), VreC, 

where C = {x £ : 3t £ (0,1) s.t. \x — te^\ < |} is a cone in -Bi(ejv). We combine the symmetry 

property and decreasing monotonicity in our proof of the nonexistence. This phenomena is due to 
the nonlocal characteristic of fractional Laplacian that requires the functions to be in L^ oc (M N ). 
Finally, our interest is to study the asymptotic behavior of u Q)P as a goes to 1“. 

Theorem 1.3 Assume that a £ (0,1), p > and u a)P is the unique weak solution of problem 
m- Then {u QiP } a vanishes as a —> 1 . 

Remark 1.2 (i) forp > a sequence of barrier functions, which converge to 0 locally in -Bi(O), 

could be constructed directly to control {u a ,p}ai 

(ii) for < p < Theorem 3.1 in )15j / is involved to control {u QtP } a ; 

(in) for p < 1 + jj, there exists a p £ (0,1) such that p < 1 + ^ for a £ (a p , 1), there is no weak 
solution for problem 11.11) from Theorem 11.21 

(iv) for p £ [1 + jt, open for the limit of {u a)P } a as a —> 1~. 

In Section 2, we treat the problem (11.141) . When the Dirac mass concentrates at point — se/v 
away from D, we build the existence, uniqueness weak solution u s of (11.141) and show how the Dirac 
mass is transformed into the nonhomogeneous term. In this case, the test functions space could be 
improved into X QjBl ( eAr ), since the solution has no singularity in D. 

In Section 3, we give a detailed account of the procedure enabling us to move the singular points 
{— sejv} to the origin. The first difficulty arises from the fact that G Q [r s ] blows up everywhere as 
s —>• 0 + [see Lemma 3.1], that is, there is no solution of 


(—A) a u = 0 in B i(e N ), 

u = 5 q in M N \Bi(e]\r). 


(1.15) 
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Therefore, we have to resort a barrier function, that is the minimal classical solution of 

(-A) a u + u p = T 0 in Bi(ejsr), 

u = 0 in Bf(ejy) \ {0}. 

In order to control the limit of {u s } near dB\, especially near the origin, some typical truncated 
functions have to be constructed carefully. The second difficulty comes from the proof of the 
uniqueness. We proceed by contradiction, assuming that there is two solutions and we will show 
their difference could be improved the test function from Cfi° into X a Bl ( ew ), this enables us to use 
Kato’s inequality US Proposition 2.4] and to conclude. 

Section 4 is devoted to blow-up case. The difficulty is to obtain the blow-up everywhere in 
-Bi(ejv) just from a lower bounds of u s , see Lemma 4.1. To overcome it, we combine the symmetric 
of the domain and resort the symmetry result of u s and then one point blowing up leads to blowing 
up every where in B i(ejv)- 

Finally, we analyse decay approximation of the weak solution for problem (11.11) when a —> . 
For p > the first challenge is to construct a sequence upper bounds that converges to zero. To 
this end, we have to study lim^^- (—where < h (J (x) = \x\~ a and then use proper parameters 
to construct the bounds. For 7 ^ 7 } < p < could not be used to construct properly the 

upper bounds, and then we use some argument of m- 


2 Dirac mass concentrated at {—scat} with s £ (0,1) 

The purpose of this section is to introduce some preliminaries. First let us state an important 
Comparison Principle. 

Theorem 2.1 | 6 l Theorem 2.3] Let u and v be super-solution and sub-solution, respectively, of 

(—A ) a u + h(u) = f in O, 

where O is an open, bounded and connected domain of class C 2 , the function f : O —>• R is 
continuous and h : M —)• K. is increasing. 

Suppose that v(x ) < u(x), Vx £ O c , u and v are continuous in O. Then 

u(x) > v(x), Mx £ O. 


Now we investigate the weak solution of 

(— A) a u + u p = 0 in Hi(ejv), 

u = 5- S6n in R N \Bi{e N ), 

where s £ (0,1). To this end, we construct a sequence of C 2 functions to approximate the Dirac 
measure. Let go : R N —> [0,1] be a radially symmetric decreasing C 2 function with the support in 
B i(0) such that f RN go(x)dx = 1. For any n £ N and s £ (0,1), we denote 

g n (x) = n N go(n(x + se N )), Mx £ R N . 


Then we certainly have that 

9n &-se N as n —> + 00 , 

in the distribution sense and for any s > 0, there exists N s > 0 such that for any n > N s , 


supp(5n) C Bs(-se N ). 

In order to investigate the solution of (12.11) . we consider the approximating solution w n of 

(— A) a u + u p = 0 in L?i(ejv), 

u = g n in R N \Bi(e N ). 


( 2 . 2 ) 
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Lemma 2.1 Assume that p > 0 and {g n } is a sequence of C 2 functions converging to 5- seN with 
supports in B^—seiy). Denote that 

g n (x):=c N: a[ - — 9n %+ 2 a d ' ! A V x £ Bi(e N ). (2.3) 

JRN \x - y\ N+2a 

Then problem h2.S\) admits a unique solution w n such that 

0 <w n <G a [g n ] in Bi(e N ). 

Moreover, the function w n := w n XB 1 {e N ) ^ Le unique solution of 

(-A ) a u + vP = g n in Bi(e N ), 

(2.4) 

u = 0 in R N \Bi(e]\r). 

Proof. The existence and uniqueness of solution to problem (12.21) refers to [JS], Theorem 2.5]. For 
n > N s , we have that supp(g n ) C £>s(—se jv) and then g n £ C 1 (Bi(ejv )) and 

W n = w n - g n in R N . 


By the definition of fractional Laplacian, it implies that 

(-A ) a w n (x) + w n (x) p = (-A ) a w n (x) - (-A ) a g n {x) + w n {x) p 

9n(z) 


CN ’ a J RN \z-x\ N +^ 


dz = g n {x). 


Therefore, w n is a classical solution of (12.41) and 

Wn < <S«[5n] in Bi(e N ), 

which implies that 

w n < G a [g n ] in B 1 (e N ). 

The proof ends. □ 

We remark that w n is the classical solution of (12.41) . then by Lemma 2.1 and Lemma 2.2 in [9], 


we have that 


J B i(ejv) 

Here (12.51) holds even for f 

Lemma 2.2 Let { g n } be defined in 
C 6 (Hi(ejv)) for 6 £ (0,1), where 

r s (a;) = 


[u; n (-A)“£ + w p n f]dx = [ fg n dx , Mf £ C™(Bi(e N )). 

J Bi{e N ) 


(2.5) 


then g n converges to T s uniformly in B i(e/v) and in 
c N,a _ tn,N 


| a; + se^\ N+2i 


Vx £ M JV \ {—scat}. 


( 2 . 6 ) 


Proof. It is obvious that g n converges to T s every point in Hi(ejv)- For x,y £ F>i(e/v) and any 
n £ N, we have that 

1 1 


\g n {x) - g n {y)\ = c N , a \ 

Je 

J B 


< CN, ( 


< 


s (_ sejv) \x - z\ N + 2a I y - z\ N + 2c 

||x — z\ N+2a — \y — z\ N+2a 
,{-se N ) \x - z\ N + 2a \y - z\ N + 2a 

|x — 2:|' V+2 " -1 + | y — z \ N + 2a ~ l 


-}g n {z)dz\ 

■g n {z)dz 


CN,a(N + 2a)\x -y\ 

J B : 


,(-«w) |x — z\ N+2a \y — z\ N+2c 


g n {z)dz 


< c 3 |x — y\ 
= c 3 |x — y\, 


Bs (— se N ) 


g n {z)dz 
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where C 3 > 0 independent of n. So {g n }n is uniformly bounded in C 0 , 1 (B i(eTv)). Combining the 
converging 

g n -t T s every point in Bi{e N ). 

We conclude that g n converges to uniformly in B\(eN) and in C e (Bi(eiy)) for 6 £ (0,1). □ 

Proposition 2.1 Assume that p > 0, s £ (0,1) and T s is given by if 2. 6 \) . Then problem 112.1\) 
admits a unique weak solution u s such that 

0 < u s (x) < G a [T s ], x £ Bi(e N ). (2.7) 

Moreover, u s := u s XB 1 (e N ) the unique classical solution of 

(-A) a u + u p = T S in -Bi(ejv), 

u = 0 in M. N \Bi(eN). 

Proof. Existence. It infers by Lemma 12.11 that the solution w n of m satisfies that 

0 <w n <G a [g n ] in Bi(e N ). (2.9) 

By Lemma 12.21 we have that g n converges to T s uniformly in -Bi(ejv) and in C e (Bi(eisr)) with 
0 £ (0,1). Therefore, there exists some constant C 4 > 0 independent of n such that 

&a[g n ](z) < 7 - C4CN \ a N+2a < C 4 C N: aS~ N ~ 2a , Vx € Bi(e N ). 

Thus, 


~N-2a 


\ w n\\mBi(e N )) < CiCN,aS N 2 a \Bi(e N )\ 


ll^n ||l oo (S 1 (ejv)) < C±CN,aS 

By [ 281 Theorem 1.2], we have that 

H^ S 'llc“(Bi(e J v)) - c 5[IW|L°c,(Bi(eAr)) + ll r s llL°°(B 1 (e JV ))] 


< C 6 [s 


—N—2a 


_|_ s -(iV+2a)pj 


( 2 . 10 ) 


for some C 5 , c§ > 0. 

In order to see the inner regularity, we denote Oi the open sets with z = 1,2,3 such that 

C>i C Oi C C >2 C 0 2 C C> 3 C 0 3 C Bi(e N ). 

By nn Lemma 3.1], for f3 £ (0, a), there exists C 7 ,cs > 0 independent of n such that 

\\ w n\\cP{ 0 2 ) ^ c 7[\\ w n\\L l (Bi(e N )) + IU°°( 03 ) + \\ w n IIl°°( 0 3 )] 

< C8[s _JV_ 2 “ + s -( JV + 2a )P], 

It follows by [28, Corollary 2.4] that there exist cg,cio > 0 such that 

\\ w n\\c 2a +P(Oi) < eg [ 11 11 (ejv-)) + 11 w n 11 CP{0 2 ) + \\ w n\\cP{0 2 )\ 

< Ciofs'^ -2 " + ( JV + 2 a )P], 

Therefore, by the Arzela-Ascoli Theorem, there exist u s £ C' 1 2 “ +e in for some e £ (0 ,/3) 

and a subsequence {w nk } such that 

w nk —> u s locally in c 2a+e as —> 00 . (2-12) 
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Passing the limit of (12.51) with £ S X aBl ( eN ) as n k oo, we have that 

[ [u s (-A) a ^ + u p s C]dx = [ £(x)r s (x). 

JB-\(eiv) JBi(eAr) 


(2.13) 


Moreover, since —>• and —>■ T s uniformly in 5i(ejv) as n —> oo, then it infers that 

0 < u s < Gaps] in -Bi(ejv)- 

Uniqueness. Let v s be a weak solution of (12.21) and then cp 3 := u s — v s is a weak solution to 

(-A) a tp s + u p s - v p s = 0 in 2?i(ejv), 

= 0 in M JV \5i(ejv)- 

By Kato’s inequality [9l Proposition 2.4], 

/ l<^l(-A)"£ + / [< - v p ]sign(u s - u s )£o?x = 0. 

JB-\(eiv) JB^ieiv) 


Taking £ = G Q [1], we have that 


/ [uf — uf]sign(u s — v s )£ t dx > 0 and / jy5 s |dx = 0 
J B\(eN) l-Bifeiv) 


Bi(en) 


then tp s = 0 a.e. in B i(ejv)- Then the uniqueness is proved. 

Furthermore, we see that w n = w n — g n is the unique classical solution of 

(-A)“u(x) + u p (x) = g n (x), Vx € Bi(e N ), 
u{x) =0, Vx € B\{eN) c 


(2.14) 


and w n converges to u s uniformly in B i(ejv)- By Stability Theorem j5j Lemma 4.5] and (12.121) . 
u sXB\( epf) is the classical solution of (12.81) . □ 


3 Proof of Theorem 11.11 

In this section, we prove Theorem II.II bv moving the points {—sejv} to the origin. To this end, we 
need derive more properties for u s , where u s is the unique weak solution of (EH). 

Lemma 3.1 Let p > 0, s £ (0,1) and u s be the unique weak solution of \2.1\) . Then the mapping 
s i —y u s is decreasing, that is, 

u si > u S2 if si < s 2 . 

Proof. By Proposition l2.il u s := u s XB\(e N ) is the unique classical solution of (I2.8|) . 

We claim that the mapping: s i->- T s is decreasing. For x £ B\{e^) and si < s 2, we observe 
that |x + siejsr | < |x + S2ejv|, then T si (x) > r s2 (x). The claim is proved. 

Therefore, for s\ < S 2 , u Sl and u S2 are super solution and solution of (12.81) replaced r s by r s2 , 
then it infers by the Comparison Principle that u si > u S2 in .Bi(ejv)- □ 

Lemma 3.2 Let s £ (0,1) and denote 


G a [T s \(x)=f G a (x,y)T s (y)dy , ^x € B 1 (e N ). 

J Si(ejv) 

Then 

lim G Q [r s ](x) = +00, Vx £ -Bi(ejv)- 

s->0+ 
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Proof. Using | 8 [ Theorem 1.2], it follows that 


G a (x , y) > min 


c n c u p a (x)p a (y) 


x,y £ Bi(e]y), 


\x — y\ N ~ 2a ’ \ x — y\ N 

where c\\ > 0 dependent of N,a. Now for x £ Bi(ejy) and y £ B i(ejv) H B\^\ (0), we have that 

4 

ciip a (x)p a (y) 


G a (x,y) > 


and 


G a [T s }(x) > 


mm 


cn 


\x - y\ N 

cnp a (x)p a (y) 


CN,a 


> 


iB^en) l\x-y | y_2a ’ \x~y\ N 

f c 1 ip a (x)p°(y) C N) g 

!B 1 (e N )nB^ i (0) \x-y\ N \y + se N \ N+2a 


| y + se n\ n+2 ° 
dy 


dy 


> \cup a (x)\x\ N [ 

5 Je 

—> Too as s —> 0 + . 


CN,aP a {y ) 

Bi( eA r)nB N (0) \y + se N \ N+2i 

~T~ 


:dy 


The proof ends. □ 

From Lemma 13.21 it is informed that the limit of G a [T s ] as s —> 0 + can’t be used as a barrier 
function to control the sequence {u s }. So we have to find new upper bound for sequence {u s }. 


Proposition 3.1 Let Tq be defined in hl.l0\) and 


2 a 

P>1 + V’ 


(3.1) 


then problem 


(-A) a u + u p = T 0 in Bi(e N ), 

u = 0 in L?f(ejv) \ { 0 }. 
admits a minimum positive solution uq, that is, uq < u for any nonnegative solution u of 13.2\) 


(3.2) 


Moreover, 


and 


lim u o(x) = 0 

aieiJi(ejv), x— >-9Bi(ejv)\{0} 


1 N+2a N+2a 

—t p < u 0 (teN) < cut p , t £ ( 0 , 1 ), 
c 12 


(3.3) 

(3.4) 


where c \2 > 1 is independent of a. 


Proof. The existence of solution to \3.2\) . It implies by Lemma 13.21 that the mapping s i—>■ u s is 
decreasing in B i(ejv), where UsXbu e N ) is the solution of (12.811 . So what we have to do is just to 
find a super solution U of (13.211 such that uq < U in i?i(ejv)- To this end, we consider the radial 
function 


®<t{x) = 


1 


Vx £ 


sN 


\{ 0 }, 


where a £ [0, N). By scaling property of dr, (also see [13]) we know that 

c(cr, a) 


(—A)“4> cr (x) = 


\a-\-2a ■ 


(3.5) 


(3.6) 
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where c(a,a ) £ M. Now we choose a = <to = then <7o £ (0, AT) if p > 1 + Therefore, 

there exist some k > 1 dependent of |c(<7o,a)| and cat iQ! but independent of n such that 

U(x) = k$ ao (x) (3.7) 

is a super solution of (13.21) . Thus, U £ L 1 (i?i(eAr)) and by Theorem 12. II we have that 

0 < u s <U for any s £ (0,1). (3-8) 

For any x £ \ {0}, uq(x) := lim s _ > .o+ u s (x) < U(x) < +oo. Following the same argument 

of Proposition 12.11 we can prove that uq is a classical solution of (13.21) . Furthermore, uq is the 
minimum solution of (13.21) . 

Proof of FOI1 . Let x £ dBi(e N ) \{0}, it'i = dBi(e N )DB^y 8 (x) and K 2 = dBi(e N )DB^ /2 (x). 
Let O be an open and C 2 set such that 

B\ (civ) n B\ s \/a{x) C O C Bi(eiv) H S^|/ 2 (x). 


Then we see that 

it i C dO and dOnK 2 = 0. 

We would like to find a super solution of (13.21) in O with vanishing boundary value in K\ for 
any n. Denote 

V x = Uri + XVo, 

where U is defined (13.71) . rj is a C 2 function such that 


if * e B \x\/^ X ^ 

if x £ B\x\/ 8 (x). 

and Vo is the solution of 

(—A ) a u = 1 in O, 
u = 0 in O c . 

Since U(1 — if) is C 2 in M N , then there exists C 13 > 0 such that 



|(—A)“[/(l — rj)\ < C 13 in 0 . 

Thus, there exists C 14 > 0 such that |(—A)“C/?/| < C 14 in O. Choosing Ao > 0 suitable, we have 
that for A > Ao, 

+ in O. 

Moreover, since rj = 1 in M. N \ O, then V\ > U > u s in M N \ O. By the Comparison Principle, we 
have that for any s £ ( 0 , 1 ), 

u s < V\ in Bi(e N ). (3.9) 

which implies (13.31) . 

Proof of d 3. Iff) . For t £ (0,1), denote that 


, . i JV+2a (X~te N 

V t (x) = c p N t p V B 


t 


x £ 


iJV 


where Vb is the solution of 


(— A) a u = 1 in -Bi(O), 

in R n \Bi( 0). 


(3.10) 


u = 0 
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It deduces by m that x E Bt (tej\f), 


— N-\-2a. r\ — AT 0 

(-A ) a V t ( x ) = c p Na t~—~ 2a ^ ~ p *~ N ~ 2a 


and 


-iv,o'' 

1 


< c nJ~ 


> 


Cl5 


\x + teN\ N+2a t N+2 "’ 

where C 15 > 0 independent of t. Then there exists some constant v E (0,1) such that 
(-A) a (zW)) + ( uV t ) p < < (-A)% + u p 0 in Bt_(te N ) 

and vVt = 0 < no in Bt(te n) c - By applying the Comparison Principle, we have that 

4 

vVt < uq in M. N , 

which implies 

1 _JV+2q - _JV+2q 

uq(x) > vclr t v min V B > c 1Q c p N \x\ p , V.x E Bt (tej\r), 

B 1(0) ’ s 

:2 


(3.11) 


for some constants ciq > 0 independent of t. We complete the proof. □ 

Now we are in the position to prove Theorem ll.il 
Proof of Theorem II. 1L By Proposition 13.11 there exists minimum solution uq of (13.21) . Using 
the maximum principle argument, we have that for s E (0,1), 


u s <u 0 in -Bi(ejv). 


Since 


u 0 E L 1 (Bi(e N )), E L 1 (Bi(e N ), p 2 dx), 
where p(x) = dist(x, ctBi(ejv)), then there exists uq < uo such that 

u s —>■ uq in .Bi(ejv) and in L l (Bi{eN)) 

and 

u p —>• Uq in Bi(cn) and in L}{B\(eN),p 2 dx). 
Passing the limit in identity of (12.131) as s —> 0 + , we have that 


'Si(ejv) 


[uo(-A) Q e + u p Qdx = [ £T 0 dx, V£ E C'o°(-B 1 (ej V ))- (3.12) 

J Bi(bn) 


By the same argument, we have that ho is C 2 locally in L>i(ejv) and uoXB 1 (e N ) is the minimum 
solution of (13.21) . Thus, we have 

u 0 = u 0 in Bi(e N ). 

Uniqueness. Let vq be a weak solution of (12.21) and then ipo := uq — vq is a weak solution of 

(-A) a </? 0 + u p Q - v% = 0 in Bi(e N ), 

<A) = 0 in R N \Bi(e N ), 

that is, 

[ ( Po(-A)°‘£dx+ [ (y% - i%)£dx = 0, £ C™(Bi(e N )). (3.13) 

JB i(ev) Is,(e w ) 
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In this definition of weak solution, we can not apply Kato’s inequality f9', Proposition 2.4] directly 
due to the stronger test functions space. So we have to improve the regularity of tpo- In fact, for 
any £ G X a Bl ( eN ), let {£&} be a sequence of nonnegative functions in Co°(-Bi(ejv)), such that 

£fc->£ and |(-A) Q £ fe | < 2|(-A)“£| in Bi{e N ). 


Passing the limit of (13.131) with £& as k —> oo, we have that (13.131) holds for any £ G X Q)Sl ( ejv ). If we 
choose a sequence {£&} C which converges to sign(rtQ — Vq) ■ G atBl (e N ) [1] and then there 

exists eg > 0 such that 


/ K 

J Bi(ejv) 


UqI G a,B 1 (e J v)[ 1 ] da; < ^611 ^011L 1 (Si (ej\r))> 


therefore, \u p Q — Vq\ £ L 1 (Bi(ej\r), p a dx). 

Now we follows the Kato’s inequality [2, Proposition 2.4] to obtain that 


/ l^o|(-A) Q £+ /" [ug - «g]sign(«o - «o)£ds = 0, V£ G X QiSl(ejv) . 


' B\ (ejv) 

Taking £ = G Q [1], we have that 


/ [?Xq — UQ]sign(no — uo)£dx > 0 and / |</?o|dx = 0 

J Bi(e N ) J Bi(e N ) 

then <^o = 0 a.e. in .Bi(ejv). Then the uniqueness is proved. □ 


4 Nonexistence 


In order to prove the nonexistence of weak solution to m for p G (0,1 + ^], we will prove that 
the solution u s of problem (12.11) blows up in Bi(ei\r) as s —» 0 + . To this end, we first define a cone 
by 

C = {x £ M. n : 3t G (0,1) s.t. \x — tejv| < -}■ (4-1) 

8 

We observe that C C Bi(ejv). 

Lemma 4.1 Assume that p > 0, s £ (0,1) and u s is the unique solution of problem k2.1\) . Then 
there exists cyj > 0 such that 

lim uJx) > cn\x\~ N , Vx G C. 

s— >o+ 


Proof. By Proposition 12.11 we have that u s := u s x Bl ( eN ) is the unique classical solution of (12. 8 j) . 
For given t £ (0,1), we recall 


V t (x) = t~ N V B 


x — tepf \ 
t 


where V B is the solution of (13.101) . 
We see that for x £ Bt(teN), 

4 


(-A ) a V t (x) = t~ N ~ 2a 


and choose s < |, then 

, , , 3 

|x + seAr| < |x| + s < —t, 

which implies that 

_ 1 _ > cis 

\x + seN\ N+2a ~ t N+2a ’ 
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where cis > 0 independent of t. Then there exists some constant v > 0 such that 


and vVt = 0 < u s in Bi(tejy) c , by applying the Comparison Principle, we have that 

4 

Vt <u s in M. N , 


which implies 

u s (x ) > vt~ N min Vb > cig\x\~ N , Vx G Bt(tei\r), (4-2) 

Bi(0) s 

12 

for some constants C 19 > 0 independent of t and s. Combining the increasing monotonicity, we 
have that 

lim u s (x) > c\Ax\~ N , Vx G C. 
s—>■ o+ 

Lemma 4.2 Assume that p > 0, s G (0,1) and u s := u s Xb\( e N ) is the minimal solution of \2.A) . 
Then 

(i) u s (x',xn) is symmetric with respect to x' and decreasing with to r' := \x'\ for any xjv G (0,2); 

(ii) u s (x',xn) is decreasing in xn for xn G (1,1 + \J 1 — |x 7 | 2 ). 

Proof. By applying the same procedure of step 1 and step 2 in proof of Theorem 1.1 in nn , we 
have that 

u s {x' ,xjv) = u s (|x 7 |,xtv) for x = (x 7 ,xjv) G -Bi(ejv) 

and u s (r,X]sr) is decreasing with r = |x 7 |. Using the same argument from the other side, we conclude 
that u s (x',xn) < u s (x'. 2 — xjv) for xjv G (1,2), that is, u s (x',Xjv) is decreasing with xjv G (1,2). 

□ 


Proof of Theorem |1.2) For p = 0, the solution u s of m satisfies 

u s = G a [r s ] — G a [l]. 

By Lemma [3721 we know that G a [-] blows every where in -Bi(ejv), but G a [l] is bounded uniformly. 
For 0 < p < 1 + our proof is divided into two steps. 

Step 1. we prove that 

lim u s (x) = +oo, x G B := {(x 7 ,xjv) G B\(e^) : xjy > 1}. 

S-S- 0 + 

If there is x = (xo,to) G B such that lim s _ >0 + u s (x) < +oo, then by Lemma 14.21 we have that 
lim u s (x) < lim u s (x) in B + := {x = ( x',xn ) G ^i(ejv) : lx 7 1 > |x 7 |, xn > to}. 

s-> 0 + s-> 0 + 

Choose a nonnegative function £o £ Cq° with support in B + , then for any x G Bi^e^) \ B + , 

i a t \ - dv f £o{y)-£o(x), 

( A) 6(*) CN'oP.Y.J^ |x_ y | N+ 2 * d y 

f Zo(y) J 

CN ’ a J B+ \x — y\ N+2a dy 

< -c N ,a\x-2e N \~ N - 2a [ £o(y)dy, 

J B+ 

therefore, there exists C 20 > 0 such that 


(—A)"^ 0 (^) < -C 20 , Vx G #i(ejv) \ B + . 
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Moreover, there exists C21 > 0 such that 

|(—A)"£o(x)| < c 2 i, Vx G B + . 
Since C C (.Bi(ejv) \ -£?+)> then by Lemma 14.11 we have that 

f u s (-A) a £ 0 dx < [ u s (-A) a ^ 0 dx + [ 

JBUefj) JC JE 


Bi(e N )\(B + UC ) 


u s (-A) a ^ 0 dx 


While 


and 


+ [ Us\(-A) a £, 0 \dx 

JB+ 

< — C20 L u s dx + C21I-B+I 

—>• —00 as s —> 0 + . 

[ v%£odx\= [ u p £ 0 dx <u p (x) max.£ 0 \B+\ 
J Si few) ./£}+ RJV 


L 


£o(z) 


Bi(ejv) l x + se ^V 


iV+2a 


dx = 


£o{x) 


>B+ 


\x + se7v|- /v+2< 


-dx < l-B+l max£0, 


where |x + seAr| > 1 for x G -B + . Then taking s > 0 small enough, we obtain a contradiction with 
the identity 


/-Bi(eiv) 

Therefore, 


/ [it s (-A) Q £ + u p i]dx = c N , a [ 
JBa (e/v) «/£ 


f(®) 


Bi(ejv) |^ + Se A r| iV + 2 ' 


lim tt s (x) = +00, Vx G £>. 

s->0+ 


-dx, 


Step 2: We claim that lim s _ >0 + u s {x) = 00, x G -Bi(ejv)- By the fact of lim s _^ 0 +(x) = 00, 
x G <8, then letting x = (i, 0 ■ • •, |), for any n > 1, there exists s n > 0 such that s n —> 0 as n —> +00 
and 

u Sn (x) > n, 

then applying Lemma 14.21 we have that 

u Sn > n in B 0 = jx = (x',x N ) G B : |x'| < ^, 1 < x N < . 

For any xo G .Bi(ejv) \ 8, there exists r\ > 0 such that B ri (x 0) C 5i(ejv) \ 8. We denote by 4> n 
the solution of 


(-A ) a u + u p = 0 

in B ri (x 0 ), 


u = 0 

in B^ 1 (xo) \ Bo, 

(4.3) 

u = n 

in £ 0 . 


e that 



u s„ — 0n 

in Si(ejv). 

(4.4) 


Let (f n = <t> n - nxB 0 , then ip n = 4> n in B ri (x 0 ) and 

-A 

L 


( — A) a ip n (x) + tp p (x) = (-A)“0„(x) - 7i(-A) a XB 0 (x) + <%(x) 

dy 


= n 


|y - x\ N+2a ' 


Vx G B ri (x 0 ), 
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that is, ip n is a solution of 


(-A ) a u + u p = n [ 

Jb 


u = 0 


dy 

B 0 1y - x\ N + 2a 


in B ri (x 0 ), 
in B^(x o). 


( 4 . 5 ) 


By direct computation, 


1 


< 


L 


dy 


< c 23 , Vx e B ri (x 0 ), 


c 23 Jb 0 I V ~ x\ N+2a 

for some C23 > 1. 

Let 77 1 be the solution of 

(—A) a u = 1 in B ri (x 0 ), 

u = 0 in B r n (x'o) 

and then (2^) 33 max 771 • 771 is sub solution of ( 14 . 51 ) for n large enough. Then it infers by Theorem 
12.11 that 

- 71 I 

Vn > (-- ) p max771 • 771, Vx <G B n (x 0 ), 

2 c 23 


which implies that 
Then by ( 14 . 41 ) . 

Since xq is arbitrary in .Bi(ejv) \ B, it implies that lim n _ ) . 00 + u Sn (x) = 00 in .Bi(ejv). □ 


71 1 

K>{~—) p max77i • 771, Vx £ B Ti {xq). 

2c 23 


lim u Sn (x 0 ) > lim </> n (x 0 ) = 00. 

s->0+ n->oo 


5 The solutions vanishes as a —> 1 


5.1 The case of p > 

Lemma 5.1 (i) Let N > 2 and cjv iQ , define in Then 


lim 
a— »1 _ 


Cat, a 
1 — a 


AN 

|5^’ 


where \S N 1 | denotes the (N — 1)—dimensional measure of the unit sphere S N 1 . 
(ii) Let N >2, then for any f € 


lim (—A) a / = -A/. 

a— 


Proof. The proofs of (i) and (ii) see Corollary 4.2 and Proposition 4.4 in m respectively. 


□ 


Proposition 5.1 Let be defined by f3.5\) with a € (0, N) and 


Then 

Moreover, 

and 


c(a,a) = -- f 

2 Jri 


\z + eisf\ a + \z-eisr\ a - 2 
\z\ N+2a 


dz. 


(-A)“<Mx) = M Virt»\{ 0 }. 


lim c(a, a) = (N — 2 — a)a 


lim (—A)“ 3 > CT (x) = — ‘f 4 _ 9 Cr ^ <7 , Vx G \ { 0 }. 


a —>T — 


| X |(T + 2 


( 5 . 1 ) 


( 5 . 2 ) 
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Proof. By direct computation, we have that 


(-A)"®„(x) = Vx£K w \{0}, 


\x 


(5.3) 


where 


z{a,ot) = ~ [ 
4 Jr 1 


i i ,_ i 

\z+e N \ a h-ejyj 
\ z \N+'2a 


- 2 


-dz. 


On the other hand, it has been proved in jl] that 

c(<7, a) = 0 if a = N — 2a. 

Now for any R > 1, let t/r : —>• [0,1] be a nonnegative function such that tjr = 1 in 

B R { 0) \ Bj_(0) and pr = 0 in B% R ( 0) U Bj_( 0), then pR^ a € Cq°(R n ). By Lemma 15.11 (ii), we 

R 2 R 

have 


lim (-A) a (p R ^ a )(e N ) = (N - 2 - a)a 

ol — y 1 


and 


{—A) a (r] R & a )(e n) = (—A) Q d> (T (eAr) -cy a f 

Jm 

A straightforward computation implies that 

(1 -Vr(z))\z 


(1 -m{z))\z\- 
| z - eN\ N+2c 


-dz. 


(5.4) 


0 '' 


< 


/ 

J R N 

I 

Jb i 


I z - ejvl^" 1 " 2 " 


dz 


— Pat V+2a 


dz 


(O)UBfj(O) I 2 “ e Af 


< C24 


/ 

JB i 


(O)UBt(O) \ z ~ e N 


\N 


dz , 


where C24 > 0 is independent of a. Thus, 


lim c/v. 

a—»1~ 


Jr 1 


(1 - 1]r(z))\z\ 
\z - ejv| v+2a 


-dz = 0. 


Therefore, passing to the limit of (15.41) as a —> 1 , we conclude that lim a ^. 1 - c(a, a) = ( N — 2 — er)cr 
and combining ([5.31) . (15.21) holds. □ 


Proposition 5.2 Let p > cr p = A+H anc ^ Uap tf ie un iq ue weak solution of hl.l\ ). Then 


u a , P < (4 1 a c N>a )p^ a in B x (e N ). 


Proof. It follows by (|5. lj) that 


c((Xp, ex ) 


(-Ar^ p (x) = g^ > vxGiym, 

where c(<7 p , q) —> N — 2 — <r p as a —> 1~. For k > 1, denote 


(5.5) 


(5.6) 


Uk — kQ ap . 

When p > we have A — 2 — a p > 0, then there exists a p £ (0,1) such that c(a p , a) > 0 for 

a € (op, 1). When p = we have a p = N — 2 < N — 2 a. Using [13], c(-,a) is C 2 and convex 
in [0,4V) and 

lim c(a, a) = — 00 
cr—yN~ 
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and for a G (0, N — 2a) 


c(a, a) > c(N — 2a, a) = 0. 

, i 

Thus, c(a p , a) > 0 by the fact of a p < N — 2a. Choosing k = (4 1 “cjv,a) p , then for any x G B i(ejv), 

( A \ a T T i rrP c N,a _ kc(<J p ,a) kP c iV,a 

l Ofc -I- jV+2a ^ \ x \cr p +2a “ l " | x |AT+2 | x |iV+2a 




|iV+2 


= 0 , 


where we used that ^ N+2 for x G Bi(ej\r). Therefore, for a G (a p , 1) we have that 

u QtP (x) < (4 1_ “cjv iQ )7<l> (x), Vx G Bi(e N ). □ 


Proof of Theorem 11.31 It follows by Proposition 15.21 that 

u a , P < (4 1- “cjv,a) p $<r p in Bi{e N ). 

On the other hand, we deduce from Lemma 15.II that 

CN,a < C25(l — a), 

where C25 > 0 is independent of a. Therefore, we have that 

1 

0 < lim u atP (x) < c 2 6 lim (1 - a)p$ rT (x) = 0, Vx G -Bi(ejv), 

a—>1~ a—>l~ 

where C26 > 0 is independent of a. This ends the proof. □ 

5.2 The case p G [jfzj, ^z§) 

We know that the unique weak solution u a , p of mn> satisfies the identity 
/ [« Q ,p(30(-A) Q £(x) + < iP (x)£(x)] dx = [ £(x)To !a (x)dx, V£ G C'^°(Bi(ejv))- (5.7) 

•liJi(ejv) JBi(ejv) 

From Lemma 15.11 part (ii), it infers that for any £ G Co°(Bi(e]y)), 

(—A) Q £ —>• — A£ uniformly in L?i(ejv) as a —>• 1~. 


Ail < 0 < A+2 _ iV+2 

AT—1 — V ^ Af-2 • °P “ p 

f/iere exists C27 > 0 independent of a such that 


Lemma 5.2 Let < p < a p = and u a:P be the unique weak solution of Then 


0 < u a>p < C27&a p in -Bi(ejv). 


(5.8) 


Proof. For < P < ]y^2> we have a p < N. It follows by (15.11) that 

(-A)«4v» = Vz 6 «" \ {0}, 

where c(cr p , a) —>• A — 2 — <7 P as a — > l - . There exists a p G (0,1) such that 

\c(a p ,a)\ < 2\N - 2 - a p \, a e (a p , 1). 

For k > 1, denote 

Uk — 
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For a £ (a p , 1) and any x £ Bi(ejv), we deduce that 

( \\ a T T i T tP c N,a ^ 2k\N — 2 — a p \ k p c iV,a 

l Uk i- U k I |jv+2a — \ x \c p +2a ^ | ;c | Ar + 2 1^1^+201 

> (feP _ 4 l ~ a c N , a - 2 N+3 ~ ap k\N - 2 - a p |) , 

which implies that there exists ko > 0 independent of a. such that for a £ ( a p , 1) we have that 

u a , p (x) < ko$ ap (x), Vx £ Bi(e N ). □ 

We have that <h CTp £ L 1 (i?i(eAr)) and 

II^Cp IImp* (Bi(0),dx) — c 28, 

where C28 > 0, p* = > 1 and M p * (L>i(0)) is the Marcinkiewicz space of exponent p*. By 

regularity results, we have that 

ll n «,pllc^ c (Bi(o)) — ° 29 ’ 

where C29 > 0 is independent of a by (15.81) . Therefore, by compactness, we only have to prove that 
0 is the only accumulation point of the sequence Let u* in L 1 (Bi(ejv)) be accumulation 

point of the sequence {u QyP } a and a subsequence, still denote {u QjP } a , converge to u*. 

On the one hand, To^ converges to 0 uniformly in any compact subset of Ri(eTv). Therefore, 
for £ € Cq^-BiIT)), we have that 


lim / £(x)ro,a(r)<Lr = 0. 

J Bi(ejv) 

Thus, we have that u* satisfies the identity 

[ [u*(x)(-A)£(x) + (u*) p (x)£(x)] dx = 0, V£ G C^iB^ex)), (5.9) 

J Bi(ejv) 

On the other hand, by regularity result, we have that u* £ C 2 (B i(eAr)) 0 L 1 (Bi^ex)) continuous 
up to boundary dB i(ejv) \ {0}, therefore, u* is a nonnegative classical solution of 


—Au + u p = 0 in Bi(ew), 

u = 0 on dBi(ei\[) \ {0}. 


(5.10) 


Then from Theorem 3.1 in m , we have that 


u* £ L°°(B i(0)) 


and Theorem 3.1 in [15], u* is a classical solution of 


—Au + u p = 0 in Bi(e]y), 
u = 0 on dB i(e N ), 

where p £ (1, (y^i). By Strong Maximum Principle, we have that 

u* = 0 in Bi(e]y). 


(5.11) 


Acknowledgements: The authors want to thank the referee for carefully reading the paper and 
making suggestions that resulted in a great improvement in the clarity of the proofs. H. Chen 
is supported by National Natural Science Foundation of China, No: 11401270 and the Project- 
sponsored by SRF for ROCS, SEM. 


18 








References 

[1] D. Balenu, K. Diethelm, E. Scalas and J. Trujillo, Fractional calculus models and numerical 
methods, Amsterdam, Word Science Publish, Nonlinear Science and complexity. 

[2] M. Bidaut-Veron, A. C. Ponce and L. Veron, Isolated boundary singularities of semilinear 
elliptic equations, arXiw.0902.0449. 

[3] M. F. Bidaut-Veron and C. Yarur, Semilinear elliptic equations and systems with measure data: 
existence and a priori estimates, Advances in Differential Equations 7(3), 257-296 (2002). 

[4] L. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, Com- 
mun. PDE 32, 1245-1260 (2007). 

[5] L. Caffarelli and L. Silvestre, Regularity theory for fully nonlinear integro-differential equa¬ 
tions, Comm. Pure Appl. Math. 62, 597-638 (2009). 

[6] H. Chen, P. Felmer and A. Quaas, Large solution to elliptic equations involving fractional 
Laplacian, Ann. Inst. H. Poincare, Analyse Non Lineaire, doi:10.1016/j.anihpc.2014.08.001. 

[7] H. Chen and H. Hajaiaj, Existence, Non-existence, Uniqueness of solutions for semilinear 
elliptic equations involving measures concentrated on boundary, arXiv:1410.267% 

[8] Z. Chen and R. Song, Estimates on Green functions and poisson kernels for symmetric stable 
process, Math. Ann. 312, 465-501 (1998). 

[9] H. Chen and L. Veron, Semilinear fractional elliptic equations involving measures, J. Differ¬ 
ential Equations 257(5), 1457-1486 (2014). 

[10] H. Chen and L. Veron, Semilinear fractional elliptic equations with gradient nonlinearity 
involving measures, J. Fund. Anal. 266(8), 5467-5492 (2014). 

[11] H. Chen and L. Veron, Weakly and strongly singular solutions of semilinear fractional elliptic 
equations, Asymptotic Analysis 88, 165-184 (2014). 

[12] E. Elshahed, A fractional calculus model in semilunar heart valvevibrations, International 
Mathematica symposium (2003). 

[13] P. Felmer and A. Quaas, Fundamental solutions and Liouville type theorems for nonlinear 
integral operators, Advances in Mathematics 226, 2712-2738 (2011). 

[14] P. Felmer and Y. Wang, Radial symmetry of positive solutions to equations involving the 
fractional laplacian, Comm. Cont. Math. 16, 1-24 (2014). 

[15] A. Gmira and L. Veron, Boundary singularities of solutions of some nonlinear elliptic equations. 
Duke Math. J. 64(2), 271-324 (1991). 

[16] Q. Guan and Z. Ma, Boundary problems for fractional Laplacians, Stochastics and Dynamics 
5, 385-424 (2005). 

[17] I. Herbst, Spectral Theory of the operator ( p 2 + m 2 ) 1 / 2 — Ze 2 /r, Commun. Math. Phys. 53, 
285-294 (1977). 

[18] E. Lieb and H. Yau, The Chandrasekhar theory of stellar collapse as the limit of quantum 
mechanics, Commun. Math. Phys. 112, 147-174 (1987). 

[19] R. Magin, Fractional calculus in bioenginering 1, 2, 3, Critical Reviews in Biomedical Engi¬ 
neering 32, 1-104, 105-193, 194-377 (2004). 


19 


[20] R. Magin, S. Boregowda, and C. Deodhar, Modeling of pulsating peripheral bioheat transf 
using fractional calculus and constructal theory, Journal of Design Nature 1, 18-33 (2007). 

[21] R. Magin, X. Feng, and D. Baleanu, Fractional calculus in NMR, Magnetic Resonance Engi¬ 
neering 34, 16-23 (2009). 

[22] R. Magin and M. Ovadia, Modeling the cardiac tissue electrode interface using fractional 
calculus, Journal of Vibration and Control 19, 1431-1442 (2009). 

[23] M. Marcus and L. Veron, The boundary trace of positive solutions of semilinear elliptic equa¬ 
tions: the subcritical case, Arch. Rat. Mech. Anal. 144 > 201-231 (1998). 

[24] M. Marcus and L. Veron, The boundary trace of positive solutions of semilinear elliptic equa¬ 
tions: the supercritical case, J. Math. Pures Appl. 77, 481-524 (1998). 

[25] M. Marcus and L. Veron, Removable singularities and boundary trace, J. Math. Pares Appl. 
80, 879-900 (2000). 

[26] M. Marcus and L. Veron, Capacitary estimates of positive solutions of semilinear elliptic 
equations with absorption, J. European Math. Soc. 6, 483-527 (2004). 

[27] E. Di Nezza, G. Palatucci and E. Valdinoci, Hitchhiker’s guide to the fractional Sobolev spaces, 
Bull. Sci. Math. 136, 521-573 (2012). 

[28] X. Ros-Oton and J. Serra, The Dirichlet problem for the fractional laplacian: regularity up to 
the boundary, ./. Math. Pures Appl. 101(3), 275-302 (2014). 

[29] V. Tarasov and G. Zaslavsky, Fractional dynamics of systems with long-range interaction. 
Commun. Nonlinear Sci. Numer. Sim. 11, 885-889 (2006). 


20 



